Lapse function for maximally sliced Brill-Lindquist initial data 



Piotr Jaranowski 
Institute of Theoretical Physics, University of Bialystok 
Lipoma 41, 15-424 Bialystok, Poland * 

Gerhard Schafer 

Theoretisch-Physikalisches Institut, Friedrich- Schiller- Universitdt 
Max-Wien-Platz 1, 01143 Jena, Germany ' 



For binary black holes the lapse function corresponding to the Brill-Lindquist initial value solution 
for uncharged black holes is given in analytic form under the maximal slicing condition. In the 
limiting case of very small ratio of mass to separation between the black holes the surface defined 
by the zero value of the lapse function coincides with the minimal surfaces around the singularities. 
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In a recent paper |1| , the authors re-obtained the Brill-Lindquist initial value solution for two electrically uncharged 
O { black holes Q using a point-mass model for the black holes followed by Hadamard's "partie finie" regularization 
, procedure in 3-dimensional space. (It can be shown that the same result appears by applying the dimensional 
regularization of Ref. ||.) In the present paper we supplement this result by the explicit calculation of the initial 
lapse function. By the time-symmetry assumption in the Brill-Lindquist initial value solution, the vanishing of the 
initial shift function automatically occurs. Therefore, all initial value metric coefficients will be known. Our result, 
| for the first time, allows the determination of the motion of test fields and particles near the initial hypersurface for 
time scales which are small against the time scale of the motion of the two black holes. Similarly, the result will 
' allow discussions of quantum gravitational problems near the initial hypersurface. Furthermore, the knowledge of the 
initial 4-metric is a step towards an (approximate) analytic treatment of the motion of two black holes starting from 
Brill-Lindquist initial value data. Finally, numerical head-on-collision simulations starting from those initial data 
may compare the numerically obtained initial lapse function with the analytic one of the present paper to check the 
accuracy of the numerical code. In our brief report, we additionally will address the question of how the surface of 
vanishing initial lapse function, which is the surface of infinite redshift for observers momentarily at rest in the initial 
hypersurface, is related in some limiting case to the minimal 2-surfaces associated with the Einstein-Rosen bridges of 
the Brill-Lindquist 3-geometry (see Ref. H ) , knowing that for the Schwarzschild black hole geometry the both types 



O 

i . of surfaces coincide (see, e.g., Ref. |4|]). 

^Jy We use units in which IGirG = c = 1, where G is the Newtonian gravitational constant and c is the velocity of light. 
. . , In the paper all vectors and their lengths are defined in the 3-dimensional Euclidean space endowed with a standard 
k*" • Euclidean metric; x is the position of an arbitrary point in this space. 

In the generalized isotropic ADM gauge [see Eqs. (7-4.22) and (7-4.23) in Ref. ||] the 3-metric takes the form 

H . / 1 \ 4 

9ij = (l + -<f>) 5 l3 + h?? (1) 

and for the canonical conjugate tt 1 ^ of the 3-metric, which in terms of the extrinsic curvature Kij of the spacelike 
hypersurfaces t = const is given by (see, e.g., Ref. Q) 

^ .</' 2 G?V m n ,J u" , ')K„„. (2) 

where g %3 is the inverse of gij and g = det(gij), the condition 

TT U EE (3) 

holds. 
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In case of time-symmetric initial conditions the black- hole linear momenta p a i (a = 1, 2) and field canonical conjugate 
ir 13 , and thus the shift function, Ni = goi, too, are zero 

Pai = 0, 7T ij '=0, Ni = 0. (4) 

Conformal flatness, i.e. 

m = + ( 5 ) 

additionally implies the vanishing of the transverse-traceless part of the 3-mctric, 

hj? = 0. (6) 

In this case, the coordinate condition (^|) coincides with the maximal slicing condition, i.e., g l] Kij — 0. Under the 
above conditions, the lapse function No (N§ — —goo + g^N^Nj — —goo, as Ni — 0) fulfills the equation 

g 1 ' 2 g i *D i D j No = ±N g l ' 2 R, (7) 

where Di means the covariant derivative with respect to gij, and R is the intrinsic curvature of the hypersurface 
t = const. The Eq. (jjj) is easily obtained from the field equation for 7r", i.e., from &7r n /dt = [see, e.g., Eq. (2.7a) 
inRef. §]. 

After substituting the ansatz 



into the Eq. (Q), we obtain the following equation for the function ip: 



1 + \<j> ) Aip = (l- ) A<p. (9) 
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To derive Eq. (|9j) we have employed the relation, valid for application on scalar functions 



and the following formula, valid when hj? = (see, e.g., Ref. ||): 

g 1/2 R = -( 1 + 

The Brill-Lindquist solution reads, 

* = *(^ + 2), (10) 

where the parameters a\ and ai are constants and where r\ and ri denote the coordinate distances between the field 
point and the black hole positions in the conformally related flat space. The parameters ot\, a,i are related to the 
bare masses mi, m,2 of the two black holes through the formulas (see Appendix A in Ref. JtJ) 



1 f n m 2 ~ m l \ 1 / m l + m 2 / "ll — m 2 X 2 



1 / mi — nii \ 1 / mi + m,2 ( TUi — x 2 



where ri2 denotes the coordinate distance of the two black holes. 
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Equation (g) is invariant against ip — > —(/>, <p — » —ip, which obviously suggests identical structures of the functions 
and ip involved in this equation. We therefore search for a solution of the function ip of the form 



ri r 2 



(12) 



where (3\ and /3 2 are some new constants. After substitution Eqs. (|T^) and (|T^) into Eq. (0) we make use of the 
regularization procedure of Ref. Q (or Ref. which relies here on replacing /(x)<5(x — x a ) by / rcg (x a )<5(x — x a ), 
where the regularized value /r C g(x a ) of the function / at its (possibly) singular point x = x„ wc define by means of the 
Hadamard's "partie finie" procedure. We expand / (x a + en) (where n is a unit vector and s > is a number) into 
a Laurent series around e = and as the regularized value of the function / at x a we take the zero-order coefficient 
of the series averaged over all directions n: 



/ + en) = ^2 fl m(n) e"\ / reg (x a ) = — j> dVL a (n 

m— — N 



The usage of the definition (|l^) boils down in our computation to simple formulas 



S(x - xi) _ S(x - xi) _ S(x - xi) 

i i u. . . - . . . 



X-Xi 



X-X 2 



Xl - x 2 



5{x - xg) _ Q <5(x - x 2 ) _ <5(x - x 2 ) 



x-x 2 



X-Xi 



xi - x 2 



(13) 



(14a) 



(14b) 



Making use of Eqs. (14) we have obtained, by comparing the coefficients at Dirac delta functions, the following 
algebraic equations for the constants (3\ and /3 2 : 



1 



a 2 
ri2 



0i 



1 



ri2 



Oil, 



(15a) 



1 



ai 
ri2 



01 



1 



01_ 
ri2 



a 2 . 



The unique solution of the both equations is easily obtained. It reads 

7*12 + Oil — OL2 



0l=OL\- 



ri2 + Oil + oi2 



(15b) 



(16a) 



02 = Ol 2 



Til + Oi2 — Oil 

m + oti + a 2 



(16b) 



In the limit ri 2 — > oo, where two completely separated Schwarzschild black holes have to arise, it holds: (3\ = ai and 
02 = OL2- If Oil = 0, then (3\ = and /3 2 = a 2 (analogously, if a 2 = 0, then (3 2 = and /3i = ai). Also, if ai = a 2 , 
then /3i = /3 2 . Thus, the solution for ■!/> does have correct limiting values. Furthermore, the post-Newtonian expansion 
of our solution for tp coincides with the corresponding part of the lapse function obtained in Ref. |(| through second 
post-Newtonian order. 

For completeness, we give the initial 4-dimensional line element. It takes the form 



ds A = - 



l-V/8 



1 + 



dt 2 + (1 + <j>/8) 4 (dx 2 + dy 2 + dz 2 ), 



(17) 



where x, y, z denote Cartesian coordinates and where the Eqs. (10), (12), and ( |16| ) [and also Eqs. (|llj) if the dependence 
in terms of bare masses is needed] have to be substituted. 

In the case of Schwarzschild black holes (there, ip = 4>), the surface N = (being there the cross-section of a 
t — const slice with the event horizon) coincides with the throat of the Einstein-Rosen bridge which is a minimal 
2-surface in 3-space. For binary black-hole systems, the surfaces Nq = need not to be identical with minimal 
2-surfaces, respectively throats, in 3-space. However, an explicit calculation has shown that the initial value surface 
No = (or, -0 = 8), coincides with the minimal surfaces around the singularities in the limiting case of very large 
separation between the black holes, i.e., we exactly have reproduced the Eq. (27) in Ref. [0] for uncharged black holes. 
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